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A mesoscopic model of amorphous plasticity is discussed in the context of depinning models. 
After embedding in a d + 1 dimensional space, where the accumulated plastic strain lives along the 
additional dimension, the gradual plastic deformation of amorphous media can be regarded as the 
motion of an elastic manifold in a disordered landscape. While the associated depinning transition 
leads to scaling properties, the quadrupolar Eshelby interactions at play in amorphous plasticity 
induce specific additional features like shear-banding and weak ergodicity break-down. The latters 
are shown to be controlled by the existence of soft modes of the elastic interaction, the consequence 
of which is discussed in the context of depinning. 


I. INTRODUCTION 

Most liquids flow as soon as they experience shear 
stress. In contrast many complex fluids (pastes, foams, 
colloidal suspensions, etc.) do not flow for shear stresses 
lower than some threshold yield limit. The rheological 
behavior of these yield-stress fluids parallels the plastic¬ 
ity of amorphous solids (oxide and metallic glasses, poly¬ 
mers, etc.). Both families of materials exhibit a rich phe¬ 
nomenology. Close to the yielding threshold, critical-like 
behaviors are observed: avalanches mm, growth of a 
correlation length scale [3], Hershell-Bulkley law [4]... In 
parallel other properties are reminiscent of glassy phe¬ 
nomena: e.g. thermal [5] and mechanical IHH] history 
dependence. In the same spirit, strain localization mm, 
a phenomenon of crucial technological interest (since it 
controls the mechanical strength) can be analyzed as an 
ergodicity break-down process: plastic activity is trapped 
in a very limited sub-region of the phase space m- 

These two phenomenological archetypes (criticality 
and glass transition) have motivated parallel modeling 
efforts. Building on trap models m designed to capture 
ergodicity breaking and aging at glass transition, Sol- 
lich et al. [131 [14] developed Soft Glassy Rheology (SGR) 
models and could associate different rheological behav¬ 
iors of complex fluids to a parameter of their model, an 
effective temperature associated to mechanical noise (see 
a recent discussion in m)- A different glassy approach 
has been pursued by Bouchbinder and Danger [16] who ex¬ 
tended the Shear-Transformation-Zone theory [TT] to ex¬ 
plicitly account for an effective temperature related to 
the slow configurational degrees of freedom of the glassy 
material under shear. 

The need to go beyond mean field description and 
understand the crucial effect of elastic interactions 
associated to the localized rearrangements (Eshelby 
events) [TTH^ responsible of amorphous plasticity has 
early led to the development of mesoscopic models ac¬ 


counting for these interactions [22]. This effort of mod¬ 
eling amorphous plasticity and/or rheology of complex 
fluids at mesoscopic scale has, since then, been very ac¬ 
tive [23H35] . As early noticed in [23], the competition at 
play in mesoscopic models between microscopic disorder 
and elastic interaction strongly reminds the physics of 
the depinning transition [36] that naturally entails crit¬ 
ical features. Recently summarized in Ref. [35], most 
features of the associated scaling phenomenology pre¬ 
dicted by depinning-like models of amorphous plasticity 
have been observed numerically [371140] and experimen¬ 
tally [HIS]. 

Noteworthily, some of the key non-ergodic features 
(e.g. aging and shear-banding [32l [4T1 02]) can be also 
recovered within the framework of the mesoscopic elasto- 
plastic models. This has raised the question of the pre¬ 
cise link with the depinning transition. In particular, the 
crucial effect of the non-positiveness of the quadrupolar 
elastic interaction induced by individual plastic events 
has been questioned. Recently Lin et al. [35] have shown 
the necessity of three independent exponents (instead of 
two for standard depinning) to account for the scaling 
properties of mesoscopic models of amorphous plasticity. 

Here we show that the specific features observed in 
elasto-plastic models are controlled by the presence of 
multiple soft modes of the quadrupolar elastic interac¬ 
tion. Note that the presence of such soft modes is not an 
artifact of lattice discretization or of a specific numerical 
implementation m- In the present perspective, shear 
bands directly result of the Eshelby interaction symme¬ 
try i.e. extended modes of plastic deformation that sat¬ 
isfy compatibility and consequently induce no internal 
stress. This property, absent in classical depinning mod¬ 
els, has dramatic effects on the stability, the dependence 
on initial conditions as well as the ergodicity properties 
of plastic yielding models. 

In the following we present in section I the details of 
the mesoscopic models of amorphous plasticity. We give 
a particular emphasis on the comparison with the models 
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of depinning an elastic manifold in a random landscape. 
The emergence of anisotropic elastic interactions associ¬ 
ated to local plastic inclusions is discussed. In section III, 
a comparison is presented between numerical results on 
strain fluctuations obtained with Mean-Field (MF) and 
“Eshelby” anisotropic elastic kernels j43l [44]. In section 
IV, a Fourier space analysis allows us to unveil the pres¬ 
ence of multiple soft modes of the Eshelby elastic interac¬ 
tions. We show in section V that this soft mode analysis 
sheds a new light on the diffusion and shear-banding be¬ 
haviors of the mesoscopic models of amorphous plasticity. 
Our main results are finally summarized in section VI. 


II. DEPINNING-LIKE MODELS FOR 
AMORPHOUS PLASTICITY 



A. A scalar mesoscopic model 

Here we restrict ourselves to a simple scalar case [24]. 
Assuming, bi-axial loading conditions, we define respec¬ 
tively for stress and strain the scalar quantities a = 
(^yy — (^xx^ ^ = ^yy “ ^xx from their tensor counterparts. 
The material is discretized on lattice at a mesoscopic 
scale ^ and is assumed to be elastically homogeneous. 
A simple plastic criterion is defined from the comparison 
between the local values of the scalar equivalent stress 
field cr(r, with a threshold 

stress cr'^[r, 5^^(r)]. The local stress a results from the 
addition of a spatially uniform external stress and 
of a spatially fluctuating internal stress due to the 
successive plastic rearrangements mediated by the elastic 
interactions. Here the local stress threshold encodes 
the disordered nature of the structure, it depends both 
on space and on the local value of the plastic strain 

From this local criterion a simple equation can be writ¬ 
ten to model the evolution of the plastic strain field: 

dts^^ (r, t)=v ((t"** + G®' *£*’* (r, t) - a" [(r, (r, t)]) 

( 1 ) 

Here the threshold dynamics is accounted for by the pos¬ 
itive part function 7^(-) such that V{x) = x if x > 0 and 
V{x) = 0 if not. 

The heterogeneity of the plastic yield stress at meso- 
copic scale is represented by the quenched variable 
The latter is defined by its average (a^) = and its cor¬ 
relations (cr^(r, z)cr^(r + 6t, f{6T)g{6eP^) 

where gives the variance. Short-range correlations are 
considered, namely, f{Sr) ^ 0 if |(5r| ^ £ and g{Sz) 0 
if \5sP^\ ^ eo- The length scale £ is given by the meso¬ 
scopic scale at which coarse-graining is performed. The 
strain scale eo corresponds to the typical plastic strain 
induced by elementary plastic events. 

Finally the internal stress is represented through 
a convolution of the plastic strain field and the elas¬ 
tic kernel associated with the reaction of the matrix 
to a unit local plastic strain: (Eshelby 

inclusion 02 ]). The properties of this long-ranged and 


FIG. 1. Sketch of a 2d amorphous material upon bi-axial load¬ 
ing. (a) The mesh is deformed according to the displacement. 
The associated strain has a reversible elastic contribution and 
an irreversible plastic contribution. The latter is represented 
according to the color scale, (b) The plastic strain field (col¬ 
ors) is represented on the un-deformed reference frame, (c) 
The plastic strain field can be represented as a d-dimensional 
manifold moving in a d -h 1 space. 


anisotropic elastic interaction are discussed in more de¬ 
tails in sub-section III Cl 

Instead of directly integrating Eq. 0 . an extremal dy¬ 
namics algorithm of the model discretized on a lattice is 
implemented [23]. Only one site (the weakest one) ex¬ 
periences plastic deformation at each iteration step. The 
external stress is adjusted accordingly. Such an algorithm 
corresponds to shearing the system at a vanishing strain 
rate and is very close in spirit to the athermal quasi-static 
protocols under conditions of imposed strain developed 
in atomistic simulations [ST] US]- 


B. From plastic yielding to depinning 

In the framework of upscaling amorphous plasticity 
from the microscopic to the macroscopic scales [46], the 
equation of evolution 0 can be understood in one of 
the two ways. First it can be seen as presented above 
i.e. as a description of the (visco-)plastic dynamics of a 
plastically heterogeneous material, discretized at scale £. 

Second, such a threshold dynamics also naturally 
emerges after coarse-graining (in the direction of motion) 
from the equation of evolution of a driven elastic manifold 
in a continuous random landscape. In order to illustrate 
this direct mapping to depinning we discuss in the fol¬ 
lowing the geometry of the equivalent manifold and the 
emergence of the threshold dynamics associated to the 
multistability of the elastic interface. 

Let us recall the equation of evolution of the over¬ 
damped motion of an elastic manifold h(x) in a random 
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landscape [36] : 

dth{r,t) = + G^’-*h{r,t) - ^ [r,h{r,t)] (2) 

Here stands for the external driving force, * h 
for the elastic restoration force and is a random po¬ 
tential such that {drU) = 0 and {drU{r,z)drU{r',z')) = 
<;‘^f{6r/i)g{6z/eo) where i and cq give the correlation 
lengths along the manifold and in the direction of prop¬ 
agation, respectively. 

The present depinning equation is very close to Eq. 
Q proposed to model amorphous plasticity. In the lat¬ 
ter the external stress plays the role of the driving force 
for the depinning, the elastic kernel associated to the 
Eshelby inclusions corresponds to the elastic restoration 
force and the disordered stress thresholds are associated 
to the random potential. 

To illustrate more clearly the direct analogy between 
deformation under shear and motion of an elastic mani¬ 
fold we give here a simple geometric interpretation. Let 
us consider the plastic strain field 5p(r) of a d-dimensional 
material. As sketched on Eig. we can define an extra 
coordinate z, orthogonal to the space variable r after 
embedding in a d + 1 dimensional space. The equation 
^ defines an elastic manifold whose propa¬ 
gation in the random landscape cr'^[r, t)] is governed 

by Eq. Q. 

An obvious difference still remains between the two 
equations. While the the depinning equation (§ models 
a continuous evolution, the equation 0 shows a discon¬ 
tinuous threshold dynamics, here encoded by the pres¬ 
ence of the p{.) function. We argue here that, far from 
being different in nature, such a threshold dynamics is a 
direct outcome of the competition between elasticity and 
disorder upon coarse-graining in the direction of propa¬ 
gation.. 

In order to give more support to the latter statement 
we resort in the following to a simple example early de¬ 
veloped in the close contexts of solid friction m and rate 
independent plasticity [48] , the over-damped dynamics of 
an isolated point driven into a one-dimensional random 
potential: 


dtx ■■ 


= -d:,W{x,y) 


(3) 


where E is a random potential such that 
{V'{x)V'{x')) = — x')/eo] where f{u) 0 

for \u\ ^ 1. Here y denotes the external driving (e.g. 
at finite velocity y = vt) and k is the strength of the 
confining potential (the stiffness of the spring driving 
the system). 

Such a system of total energy W{x,y) = k{x — y)‘^l2-\- 
V{x) is known to exhibit multistability when disorder 
overcomes elasticity. Namely if /c^/eo > 1, for ev¬ 
ery y position, one and only one position x*(^) satis¬ 
fies equilibrium and stability conditions: dxW{x,y) = 0 



FIG. 2. Phenomenology of the motion of a particle of position 
X driven by a spring of position y in a one-dimensional ran¬ 
dom potential V{x): (a) Graphical representation of the mul¬ 
tiplicity of the solutions of the equilibrium equation V'{x) = 
—k{x — y) for two positions of the spring, (b) Parametric 
represntation of the complex effective potential We // (y) seen 
by the spring and representaion (in red) of on eparticular tra¬ 
jectory. (c) Associated represntation of the force landscape. 
Jumps in the potential Weff are associated to force thresh¬ 
olds. 














































4 


and d‘^^W{x^y) > 0. An effective potential VFeff(^) = 
VF[x*(^)] can then be defined unambiguously. 

Conversely, as illustrated in Fig. that shows graph¬ 
ical solutions of the equilibrium equation —k{x — y) = 
V'{x), for k<;/eo < 1 the potential W is characterized by 
a large number of local minima and several stable posi¬ 
tions x*(^) of local equilibrium can be found for a given 
position y. 

Still, it is possible in this multistability case to resort to 
a parametric representation and to build an effective po¬ 
tential Wef^[y*{x)] = bF[x,^*(x)] associated to the mul¬ 
tiple minima. As shown in Fig. |^, the stable branches 
of this effective potential consist of series of truncated 
parabolas. Upon driving, the system jumps from one 
local minimum i to another one j as soon as a force ex¬ 
ceeds the threshold value fi = —V{xf^) associated to the 
upper bound xf^ of the basin of attraction of the min¬ 
imum i (the intersection with the next parabola). One 
obviously recovers here the phenomenology of the insta¬ 
bility inducing local rearrangements at the atomic scale 
in amorphous plasticity m- 

An example of such an (history-dependent) trajectory 
made of series of micro-instabilities is shown in Fig. 
an d|J;. A threshold dynamics thus directly emerges from 
this simple case of an isolated defect. In particular, as 
shown in Fig. it is clear that upon coarse-graining at 
scale the dynamics of jumps between basins is entirely 
controled by the series of threshold forces fi. 

The phenomenology remains unchanged when dealing 
with more complex objects like manifolds. Rather than 
the stiffness of an external device, the disorder has in 
this case to be compared with the internal elasticity of 
the manifold. See e.g. Ref. [49] for a recent discussion 
in the context of crack front propagation. Note that the 
non-regularity of the effective potential induced by mul¬ 
tistability is likely to be related to the emergence of a 
cusp in the correlator of depinning forces observed under 
renormalization m- 

The present model of amorphous plasticity this ap¬ 
pears to belong to the wider class of depinning models. 
We discuss in the next section to what extent the peculiar 
nature of the Eshelby elastic interaction associated with 
plasticity does affect the phenomenology of depinning. 


C. A peculiar elastic interaction 

The occurrence of a plastic local rearrangement in the 
amorphous structure inuces internal stresses due to the 
reaction of the elastic surrounding matrix. This results in 
a stress relaxation of the region that rearranged and in an 
anisotropic long-ranged stress field in the outer matrix. 
This elastic interaction is very peculiar. In particular, it 
is non strictly positive: the sign depends on the direction. 
The elastic interaction thus either favors or unfavors the 
occurence of future plastic events depending on their po¬ 
sition. 

The exact internal stress field obviously depends on the 


details of the rearrangement of the amorphous structure. 
A classical approximation consists in resorting to a con¬ 
tinuum mechanics analysis and in using the solution of 
the stress induced by a plastic inclusion early proposed 
by Eshelby [43]. More precisely, independently on the 
precise shape of the inclusion, only the dominant contri¬ 
bution of the internal stress in the far-field is considered. 

In the plane strain geometry considered in Eig. |^, a 
pure shear plastic inclusion induces a long-range inter¬ 
nal stress characterized by a quadrupolar symmetry. In 
an infinite medium, the dominant term in the far-held 
and the mean stress drop in the inclusion can be written 
respectively: 

GqW = , Gq(0) = • (4) 

where /i* is an effective elastic modulus, a the mean ra¬ 
dius of the inclusion and Cp the mean plastic strain expe¬ 
rienced by the inclusion. Here the subscript Q refers to 
the quadrupolar symmetry. Note that the amplitude of 
this quadrupolar elastic interaction is controlled by the 
product of the “volume” of the inclusion by the mean 
plastic strain. 

Eor the numerical implementation, bi-periodic bound¬ 
ary conditions are considered and following Ref. [24], a 
quadrupolar lattice Green function is dehned from the 
following expression in the Eourier space: 

G^(p, q) = -A [cos(40pg) + 1 ] , G^( 0 , 0 ) = 0 . (5) 

where Opg is the polar angle and (p, q) the wavevector in 
Eourier space. While the hrst term directly stems from 
the quadrupolar symmetry of the Eshelby far field ex¬ 
pression 0 , the null value of the zero frequency term 
Gq( 0 , 0 ) is required by a stationarity condition : a spa¬ 
tially uniform plastic strain induces no internal stress. 
In other words, no plastic incompatibilities are gener¬ 
ated because of the assumption of uniform elastic mod¬ 
uli. When translated to discrete form, it means that 
j ^Q(bi) = O 5 henceforth this condition directly im¬ 
poses the value of the latice Green function at the origin, 
i.e. the stress drop: 

Gq(0,0) = - ^ GQ(i,j) (6) 

= E [cos(40p,) + 1 ] 

(p,g)7^(0,0) 

The prefactor A has the dimension of an elastic modulus. 
Here it is chosen so that the local stress relaxation in the 
site that experienced a unit plastic deformation is unity: 
Gq( 0 , 0 ) = - 1 . 

In the plane shear strain geometry (invariant along the 
2 :-coordinate) illustrated in Eig. the quadrupolar elas¬ 
tic interaction Gq is positive in the directions at ±7r/4 
and negative in the directions at 0 and 7r/2. The as¬ 
sociated plastic strain field is thus orientated along the 
diagonals of the x, y plane. 
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FIG. 3. Sketch of two kinds of shear geometry: (a) Plane 
shear geometry; (b) Antiplane shear geometry. In both cases 
the strain field is invariant along the z-axis. The color scale 
gives the amplitude in the xy plane of the plastic strain fields 
(a): = Sy^y — and (b) Sy^z- A quadrupolar symmetry is 

observed in the plane shear case (a) while a dipolar symmetry 
is observed in antiplane shear case (b). 



For the sake of completeness, we also illustrate in 
Fig-i) another loading geometry: the antiplane shear 
geometry. Here the strain field is again invariant along 
the 2 :-axis but the system is sheared along the yz direc¬ 
tion so that only Uz the ^-component of the displacement 
field is non zero and the strain component of interest is 
Syz{x^y) = dyUz{x^y). Within this antiplane geometry 
early studied in Ref. [23], the elastic interaction asso¬ 
ciated to a plastic inclusion obeys a dipolar geometry: 
Gd{^) = Acos{20)/r‘^ so that the plastic strain field is 
orientated along the x direction. The specificity of this 
loading geometry will be further discussed in section |VT1 
Due to their long range character, it may be tempt¬ 
ing to approximate the “Eshelby” elastic interaction by 
a simple Mean-Field (MF) interaction [28]: GMri'^ij) = 
1/(A^^ — 1) if \r\ 7 ^ 0 and Gmf(O) = —1. The latter 
will be used (all other parameters being kept constant) 
to illustrate the expected behavior of a standard refer¬ 
ence depinning model. In the following, we compare the 
respective effects of Mean-Field and quadrupolar inter¬ 
actions on some specific properties of amorphous plastic¬ 
ity, i.e. strain diffusion and localization. In order to to 
investigate the origin of the specific effects of the “Es- 
helby” elastic kernel, we also define a weighted average 
of two propagators: Ga = (1 — a)GQ + aGMF where the 
parameter a gives the relative weight of the mean field. 
Eor moderate values of a, the quadrupolar symmetry is 
mainly preserved in the sense that the Green function 
remains strictly negative in the 0 and it/2 directions. 

III. MEAN-FIELD DEPINNING VS 
PLASTICITY 

A. Family-Vicsek scaling vs Diffusion 

We first discuss th e beha vior of the variance of the 
plastic strain where we defined the spa¬ 

tial fluctuation of the plastic strain field 
Here X and {X) denote the spatial average and the en- 


FIG. 4. Plastic train variance a nd elastic stress vari¬ 
ance vs cumulated plastic strain for a quadrupolar 
propagator Gq and a Mean-Field Gmf- A linear behavior is 
represented for comparison. 



FIG. 5. Size-dependent behavior of the plastic strain diffu- 
sivity Dq = Wq/stp^ obtained with the quadrupolar kernel 
Gq for sizes N = 32, 64, 128 with M = 1000, 100, 30 re¬ 
alizations respectively. The larger the system, the longer the 
anomalous sub-diffusive behavior. 


semble average of the variable X, respectively. We show 
in Fig. (a) the evolution of the variance with re¬ 
spect to the mean plastic strain In the context of 
depinning, as illustrated in Fig is nothing but 

the width of the propagating interface. Moreover, in the 
framework of extremal dynamics used here, the mean 
plastic strain sp^ defines a fictive time directly associated 
to the total number of iterations. It is thus legitimate 
to discuss our results in the framework of the classical 
Family-Vicsek scaling [5TH53] for interface growth. The 
latter predicts first for the interface width VF, a power 
law growth W oc up to a time scale r oc such 
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FIG. 6. Strain variance (equivalently interface width) vs cu¬ 
mulated plastic strain ^ for 8 different propagators: Gq, Ga 
(a = ± 10 -^± 10 -^± 10 -^), Gmf. 


that the correlation length ^ has reached the system size 
^(r) ~ L and beyond which saturation is obtained. 

Our numerical results are shown in Fig. for Mean- 
Field and quadrupolar elastic interactions. As expected, 
the classical Family-Vicsek scaling is recovered for the 
width obtained in the case of the Mean-Field 

depinning. In the amorphous plasticity case, the first 
power-law growth regime is recovered but, past ^ ^ L, 
the interface width Wq shows no saturation but rather a 
diffusive trend [24] . The evolution of the variance of 
the elastic stress field is also shown in the two cases. 
Here saturation is recovered in plasticity as well as in 
MF depinning. Note that the elastic stress field can be 
directly obtained from the plastic strain field from a sim¬ 
ple convolution with the propagator: = G ^ The 

fact that the diffusive trend at play with the strain field 
does not show in the stress fluctuations is a first indica¬ 
tion that strain fluctuations are controlled by soft modes 
of the elastic interaction. 

In order to characterize in more details the diffusive- 
like behavior of the plastic strain field obtained with 
the quadrupolar elastic interaction, we show in Fig.j^b) 
the evolution of the associated effective difusivity Dq = 
WqI^. This ratio is expected to be constant for stan¬ 
dard diffusion. At very short times, a plateau is ob¬ 
served; In this very early regime, plastic activity is not 
correlated yet. Then the diffusivity Dq shows a power- 
law growth. This simply derives from the fact that in 
this regime the growth exponent a is larger than unity: 
Dq = W^/^i (X 

The evolution of the diffusivity then shows a strong 
size-dependence. For small system size, a simple plateau 
is obtained, the diffusivity saturates to a constant value. 
However for larger system sizes a long decreasing tran¬ 
sient is observed before a stationary value is obtained. 
The larger the system, the longer the transient sub- 


diffusive regime. 


B. Shear-banding and plastic aging 


The nature of the elastic interaction thus strongly af¬ 
fects the evolution of the spatial fluctuations of the plas¬ 
tic strain field and in particular the existence of a dif¬ 
fusive regime. In order to get more insight on the re¬ 
spective effects of the Mean-Field and the quadrupolar 
kernels, we now show results obtained with the mixed 
kernel Ga = (1 — a)GQ + aGMF- 

In Fig. (top) the evolution of the interface width is 
shown for different (small) values of a. It turns out that 
even the lowest positive MF contribution is enough to 
recover saturation at long times. A transient diffusive 
regime appears when a tends to zero, and the level of the 
final plateau increases accordingly. But when the inter¬ 
face gets too distorted, if a > 0 the (low) MF restoring 
force eventually stops the divergence of the strain fluctu¬ 
ations. 

A negative MF contribution a < 0 has the opposite ef¬ 
fect: after a transient diffusive regime, the plastic strain 
becomes unstable and its variance diverges very fast. The 
diffusive regime thus appears to be a specific feature of 
the quadrupolar kernel. It lives on the verge of stabil¬ 
ity and any mean-field contribution to the elastic kernel 
sends the system either toward saturation or ballistic evo¬ 
lution depending on the sign of a. 

The strong effect of the MF contribution is also man¬ 
ifest in the spatial distribution of the plastic strain 
field. In Fig. (bottom) maps of the plastic strain 
are shown for a cumulated plastic strain ~ 5 for 

a = —0.01, 0, 0.01 using the same color scale. The plas¬ 
tic case (a = 0.0) shows a superposition of patterns local¬ 
ized at ±7r/4 following the symmetry of the quadrupolar 
kernel. Similar patterns survive with a positive MF con¬ 
tribution (a = 0.01) but get very attenuated (the inter¬ 
face width is much lower). A negative MF contribution 
induces conversely a strong localization behavior: plastic 
activity is restricted along a unique very thin shear band. 

As above mentioned, shear-banding can be analyzed 
as a kind of ergodicity breakdown: plastic deformation 
only visits a sub-part of the phase space m- It is thus 
tempting to analyze the present model of plastic yielding 
along these lines. In Fig. we show two-point correla¬ 
tion functions computed after various “waiting times” ^ 
(here the cumulated plastic strains): 


G{ew,£) 


{sP^{x,y,s)sP^{x,y,s^)) 

(^) V 2 ( 4 ) l /2 


(7) 


For the bare plasticity model, a striking mechanical 
history effect is observed: the larger the waiting time, 
the larger the decorrelation time. Again, the addition of 
a very small MF contribution is enough to destroy this 
mechanical history dependence. Such results are remi¬ 
niscent of recent studies of depinning lines [54] that re¬ 
vealed aging properties but only in the roughness growing 
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FIG. 7. Maps of plastic strain field obtained for a mere quadrupolar elastic interaction (b), and with a positive (a) and a 
negative (c) MF contribution a — ±10“^ for ^ 5, past the transient regime. The same color scale has been used in the 

three cases. 


stage. Here the saturation of the interface roughness is 
postponed at infinity and aging can persist forever. This 
regime is thus naturally associated to the divergence of 
the interface width. 

Note that such an aging behavior may also be observed 
in a simple diffusion process. The diffusion regime at 
play in amorphous plasticity is however highly non triv¬ 
ial [24l [39]. In particular, as shown in Fig. m (b), for 
large systems, a very long sub-diffusive transient regime 
is obtained i.e. we get Wq oc with /3 < 1 over a wide 
range of strain. This observation again supports weak 
ergodicity breakdown. The latter behavior is indeed as¬ 
sociated to sub-diffusion [55] . 


IV. FOURIER SPACE AND SOFT MODES OF 
THE ELASTIC INTERACTION 

The introduction of yet a tiny MF component has thus 
dramatic consequences on the localization behavior, a key 
feature of amorphous media plasticity. In the following, 
a rewriting in Fourier space allows one to emphasize the 
crucial role of the soft modes of the propagator in this 
phenomenon and their connection to plastic shear-bands. 


A. A Fourier representation of depinning 

In the model presented above, the “Eshelby” 
quadrupolar interaction was defined through its Fourier 
transform in order to handle periodic boundary condi¬ 
tions [24] : 

Gpq = A {— cos{4:0pq) — 1 ) = —2A ^^2 + ^ 

where Opq is the polar angle and (p, q) the wavevector in 
Fourier space. A is a constant chosen so that G(0, 0) = 


— 1. The Fourier transform of the plastic strain field is 
defined as: 


pi 

^mn 


1 


N/2-1 N/2-1 _ 

E E 4 


P=-NI2 q=-N/2 


_• 27rmp 

e ^ e 


(9) 


The Fourier components of the quadrupolar elastic inter¬ 
action is thus: 
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( 10 ) 


• 27rmp • 27mq 

Denoting Bpq = the (p, g) Fourier mode, 

we get G^^epq = XpqBpq with Xpq = -2A[{p‘^ + 
g^)]^. In other terms, the eigenmodes of the Green op¬ 
erator are precisely the Fourier modes, and the associ¬ 
ated eigenvalues are the above written Xpq. This prop¬ 
erty stems from the translation invariance of the elastic 
propagator. 

The same property also holds for the MF propagator: 

= -SmSn + (1 - SM /" 1 ) 


where N is the linear size of the system. 

Let us now discuss the eigenvalue spectrum of the 
quadrupolar interaction. One first recognizes the trans¬ 
lation mode of zero eigenvalue Aqo =0. In the classical 
depinning case (say MF, Laplacian or power-law in dis¬ 
tance) this mode is the only one characterized by a zero 
eigenvalue. It is the signature of the invariance of the 
model with respect to a uniform translation of the man¬ 
ifold along its propagation direction. 

In the quadrupolar case, a set of non-trivial eigen¬ 
modes are also characterized by a null eigenvalue. 

• 27 rp(m + n) • 27 rp(m — n) 

Namely Bp^p = ^ and Gp^-p = n with 

p G [—N/2^ N/2 — 1]\{0}. Thus there is one trivial zero 
translation eigenmode and 2{N — 1) non-trivial ones. 
















FIG. 8. Two-point correlation Ci^^e) of the plastic strain field as a function of the cumulated mean plastic strain e for 
4 “waiting times” ^ = 1, 4, 16, 64 and for 4 different propagators: Gq (a), Ga with a Mean-Field weight a = 10“^ (b), 
a = 10“^ (c) and a = 10“^ (d). A clear aging effect shows for the quadrupolar propagator Gq: the longer the waiting time, 
the slower the decorrelation. The small MF contributions in propagators Ga gradually kills the aging behavior. 


Let us rewrite the the plastic strain field in the Fourier 
basis using the more condensed form: 

= '^Cpgep^g , where Cpg = ^e^g (12) 

In order to follow the evolution of the different modes we 

now rewrite in Fourier space the argument of the 7^(-) 
function in the equation of evolution 0 : 

—_ (13) 

a^^*5p6g + G%el\ - . 

Ignoring for the moment the effect of the function 7^(-) 
in Eq. 0 we thus get by Fourier transform the evolution 
of the contribution of the different modes: 

r)c - - 

+ Xpqdcpg - [c7<=(r, (r))]p^ (14) 

This rewriting thus enables us a better understanding 
of the diffusive-like behavior observed at long times for 


the plastic strain. In real space, the spatial coupling is 
induced by the non-local elastic interaction kernel, 
while the noise term is local. In the space of eigenmodes, 
the opposite character is observed, namely the restoring 
force is local, but noise is not. Since all eigenvalues are 
null or negatives (otherwise the dynamics would be un¬ 
stable) a competition emerges between the relaxation of 
the eigenmodes induced by the elastic contribution and a 
random forcing due to the quenched noise contribution. 
In particular, at long times, the contribution of the soft 
modes becomes dominant since they are not submitted 
to relaxation. The diffusive-like behavior thus directly 
emerges from a competition between the different soft 
modes controlled by the quenched disorder. 

The strong effect of a small MF contribution to the 
quadrupolar propagator can now be re-read as the con¬ 
sequence of the opening of a gap in the spectrum of 
eigenvalues, in other words to the vanishing of the soft 
modes. In Fig. the spectra of eigenvalues of the stress 
redistribution kernel show the gradual gap opening due 
to the introduction of a MF contribution to the elastic 
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that the two systems of shear bands are not indepen¬ 
dent since the scalar products dk.d'i may be non zero. 
Here the amplitudes and roughly correspond to the 
mean plastic strain along the shear-bands dk and d'k re¬ 
spectively. In the same spirit as above, accounting for the 
non-orthogobality between the two slip systems, it is pos¬ 
sible to write the equation of evolution of the amplitudes 
of the shear-bands: 



l^<T'=[(r,£^'(r)]( 16 ) 

redfc 

red'fc 


FIG. 9. Spectrum of eigenvalues of elastic propagators: 
Mean-Field (MF), quadrupolar interaction and MF-weighted 
quadrupolar interactions. Eigenvalues are here simply ranked 
in the decreasing order. The introduction of a fraction a of 
MF opens a gap between the translational mode having a null 
eigenvalue and the other modes A < 0. The evolution of the 
gap is zoomed in the inset. 

quadrupolar interaction. The associated restoring elastic 
force brings back the model to the standard depinning 
phenomenology. 

Note that this interpretation only holds if we ignore 
the 7^(-) function that intervenes in Eq. When a long 
integration time in considered, the loading contributes to 
a positive average that allows for such an interpretation. 
However, at short time scales, the positive part func¬ 
tion unfortunately cannot be simply expressed in Fourier 
space. A similar situation appears in classical depinning 
models. The point is that for the latter ones, a long time 
integration gives a finite restoring force to any wavelength 
of manifold fluctuations. 


B. Prom Soft modes to shear-bands 


where in the present case of bands at ±7r/4, //c^/2 = 
{k — 1) (mod 2). 

As already discussed above, in absence of elastic restor¬ 
ing force in the equation of evolution, we expect the strain 
field to be asymptotically dominated by the sole superim¬ 
position of soft modes, which we interpret here as shear 
bands at ±7r/4. 

Here we obtain for the dynamics of the bands an advec- 
tion contribution due to the differnce — a^) between 
the driving force and the spatial average of the threshold 
on the whole lattice. In addition, the average along the 
bands of the fluctuating part of the thresholds and the 
inter-bands coupling introduce randomness and lead to 
diffusion. 

Note that another important souce of interactions be- 
twen bands has been neglected here. Although shear 
bands are expected to be dominant at long times, the 
short time synamics remains local. A natural conse¬ 
quence of the interplay between a local threshold dynam¬ 
ics and the non-local effects of the elastic interaction is 
the persistence of fluctuations along the bands. The con¬ 
volution of the latter with the elastic kernel is respon¬ 
sible for a mechanical noise contribution in the dynam¬ 
ics [ISIIMIET]. 


In the present context of amorphous plasticity an ap¬ 
pealing alternative representation of the soft modes is 
given by the unit shear-bands orientated along ±7r/4. 
One defines d^ such that d/c(m,n) = Sm-n-k and 
d'/c(m, n) = 5m-\-n-k where k G [—N/2 + 1, N/2 — 1] and 
6n is the Kronecker symbol. Plastic shear bands thus 
directly appear as soft modes of the quadrupolar elas¬ 
tic interaction, because of the null eigenvalue, they don’t 
induce any internal stress. 

We use this decomposition to rewrite the plastic strain 
field as: 

= E! + E] + E] 

\p\7^\q\ k k 

where the first sum gathers all modes of non-zero eigen¬ 
values whereas the two last sums correspond to combi¬ 
nations of shear bands oriented at ±7r/4. Note however 


V. FLUCTUATIONS AND AGE STATISTICS 
ALONG SHEAR-BANDS 

The interpretation of the plastic shear-bands as soft 
modes of the elastic interaction encourages us to re¬ 
examine our results from this new perspective. In partic¬ 
ular, we expect that at long times, plastic activity con¬ 
centrates along weakly interacting shear bands. A natu¬ 
ral question thus arises about the respective contribution 
of intra shear-bands and inter shear-bands fluctuations to 
the diffusive regime. This question is reminiscent of ear¬ 
lier studies showing anisotropic correlations in the plastic 
strain field [smsi]- In the same spirit, we suggested that 
the long sub-diffusive regime observed in the numerical 
results reflects an aging-like behavior. This motivates 
us to characterize age statistics inside and outside shear- 
bands. 
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We first define the mean variance of the plastic strain 
field inside the shear bands as: 

1 ^ 

^ k=l 

r 1 ^ 

redfc L rGdfc 

In the quadrupolar geometry associated to plane shear 
plasticity, the shear bands and d'k are oriented along 
the ±7r/4 directions and receive a positive stress contri¬ 
bution whenever one of their site experiences plasticity, 
hence the superscript + in the notation of the variance 
Wq . In a similar spirit we can can characterize the fluc¬ 
tuations of the plastic strain held along the directions at 
angles 0 and 7r/2 that receive a negative stress contri¬ 
bution when one of their site experiences plasticity. We 
denote Wq the variance of the plastic strain held along 
such negative stress directions. 

We show in Fig.[T^ the evolution of the global variance 
Wq of the plastic strain held as well as the variances Wq 
inside the shear-bands and Wq outside the shear-bands. 
We observe that the variance Wq of intra-shear-bands 
fluctuations are significantly lower than the global vari¬ 
ance Wq in the diffusion regime. Conversely, the variance 
Wq measured in the negative stress directions is indistin¬ 
guishable from the global variance. The inset shows the 
same data after rescaling by the mean plastic strain i.e. 
the effective diffusivities Dq = Wq /^, Dq = Wq 
and Dq = Wq /^. Here we see that the effective dif- 
fusivity within the shear-bands Dq is about two times 
smaller than the global diffusivity Dq. 

Beyond the spatial fluctuations, we can also charac¬ 
terize the temporal fluctuations. In order to do so, we 
define the local age variable Aq = nAlN‘^ that counts the 
number of plastic events ua that occurred in the system 
since the last time the site has experienced plasticity. In 
case of an homogeneous deformation, all N‘^ sites would 
be expected to experience plastic events at the same fre¬ 
quency, hence the rescaling factor 1/N‘^. It is easy to 
extend this definition to a shear-band: A'q = Ua+IN. 
Here is the number of plastic events since the last 
time a site of the band has experienced plasticity and the 
rescaling factor stems from the number N of shear bands. 
The age Aq of bands in the negative stress directions is 
defined in the very same way. 

We show in Fig. [Tq| 3 the distributions of ages 
P(log^Q), P(log^g) and P(log^g) measured in the 
diffusive regime. The age distribution of sites P(log^Q) 
peaks around unity and shows a cut-off around ten. This 
suggests that on average, the plastic activity is only mod¬ 
erately heterogeneous. 

As for the spatial fluctuations we observe that the age 
statistics of bands P(^g) measured in negative stress di¬ 
rections (outside shear-bands) is close to the global age 
statistics P{Aq) measured on individual sites. In con¬ 
trast, the distribution P(^g) of ages of the shear-bands 



FIG. 10. Quadripolar kernel, (a): variance of the plastic 
strain field inside and outside shear bands, (b): Age distri¬ 
bution inside and outside shear bands 


is shifted to larger values. A natural interpretation is that 
due to the positive stress redistribution, plastic activity 
remains trapped for longer periods within a shear-band 
(while the age of the other bands keeps increasing) before 
jumping to another one. We note in particular that the 
cut-off of the shear-band age distribution roughly corre¬ 
sponds to the duration of the subdiffusive regime. 

The spatio-temporal fluctuations of the plastic activ¬ 
ity within the shear-bands is thus clearly distinguishable 
from the backround. Still, this diference is not dramatic. 
Although the diffusivity is decreased and the duration 
of plastic activity is increased along the shear bands, 
the qualitative picture remains unchanged. Shear bands 
can survive 5-10 times longer than bands in the negative 
stress directions but the age statistics ends up converging 
toward a stationary distribution. This is for instance at 
contrast with the clear ergodicity breaking identified in 
Ref. pT] . 


VI. PLANE VS ANTIPLANE SHEAR IN 
AMORPHOUS PLASTICITY 

A potential reason for the system to escape aging actu¬ 
ally stems from the quadrupolar geometry of the elastic 
interaction at play in the present model. Since, after a 
plastic event, the elastic stress is positive along the two 
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directions at ±7r/4, it is possible to trigger another plastic 
event in a direction at 0 or tt with a sequence of two suc¬ 
cessive events at +7r/4 then —7r/4 (or the reverse). Such 
sequences thus restore some interaction between positive 
and negative stress directions. 

In this section we follow this geometric idea by fo¬ 
cussing on the case of antiplane shear geometry early 
studied in Ref. [23]. As mentioned above, in this an¬ 
tiplane geometry (defined in Fig[^), a plastic inclusion 
induces a dipolar interaction: 

O°(0) = -l (18) 

Here the soft modes are shear-bands oriented at 6> = 0 
and the negative stress directions are oriented at ^ = tt. 
In contrast with the previous quadrupolar case, no di¬ 
rect cross-talk mechanism is possible between the differ¬ 
ent shear-bands. This means in particular that if we now 
rewrite the plastic strain field as: 

gPt = CpqBp^g + , ( 20 ) 

blT^O k 

where the N horizontal bands d^k are the soft modes of 
the dipolar kernel G^, we now obtain for the equation of 
evolution of the band amplitudes: 

1 _ 

V S (21) 

red^k 

We thus get in the long term dynamics a set of of bands 
that can grow independently of each other. Again, this 
statement has to be softened to account for the effective 
noise induced by the sort term local thershold dynamics 
that restore weak coupling between the bands. 

In analogy with the previous section we show in 
Fig. [TTJi the evolution upon deformation of the vari¬ 
ances and of the plastic strain field obtained 
along the positive and negative stress directions, respec¬ 
tively, in comparison with the global variance Wd- As in 
the quadrupolar case, the variance in the negative 
stress directions is almost the same as the global variance 
Wd- The result is strikingly different in the direction of 
shear-bands. After the power-law transient, instead of a 
diffusive regime, the variance shows indeed a clear 
saturation. Along the direction of the shear-bands, we 
thus recover the classical Family-Vicsek phenomenology 
of depinning. Note however that saturation is reached 
at a much later stage ~ 10 than in the reference 
Mean-Field case ^ 0.5 (see Fig. [^. If one refers 
to the results obtained with the composite kernels Ga 
(see Fig. [^, this would correspond to small Mean-Field 
weight a ^ 0.005. 

In Fig. we show the distribution of ages in the 
antiplane shear geometry. Again, the age distribution 
(logM^) of bands in the negative stress direction is very 



FIG. II. Aging and diffusive bhaviors obtained with a Dipo¬ 
lar kernel, (a): variance Wd of the plastic strain field and 
variances inside and outside shear bands, (b): Age 
distributions P(\ogAD) of the sites and and age distributions 
P(logMj), P(logMb) bands in the positive (shear- 

bands) and negative stress directions, respectively. 


close to the age distribution P{logAD) of the individual 
sites. The case of the shear-bands is strikingly different. 
Here the age distribution P(log M^) is much older (about 
two orders of magnitude) than the two other ones. One 
recovers the same aging-like effect as for the shear-bands 
in the quadrupolar case but with a much higher ampli¬ 
tude. 


VII. CONCLUSION 

Depinning models rely on the interplay between disor¬ 
der and elasticity. While the yielding transition may be 
discussed within the framework of depinning, it appears 
that some specific properties of amorphous plasticity (dif¬ 
fusion, shear-banding) are controlled by the peculiar form 
of the quadrupolar elastic interaction. In order such fea¬ 
tures to be recovered in the framework of discrete lattice 
models, the discretized implementation of the Eshelby 
kernel has to preserve a key property of continuum plas¬ 
ticity: a unit plastic strain along any band in a direction 
of maximum shear stress (here ±7r/4) induces no residual 
stress. 
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The interpretation of the shear-bands as soft modes of 
the Eshelby elastic interaction may clarify the long de¬ 
bate about the relative importance of localized rearrange¬ 
ments and large scale shear-band like events mm 
as microscopic mechanisms of amorphous plasticity and 
complex rheology. It appears in particular that localized 
plastic events is the rule at short time scales but on a 
larger time horizon, mostly shear bands account for the 
kinematics, as these are the only displacement fields that 
prevent large shear stress build-up. 

While the present study has been concerned with the 
modeling of amorphous media plasticity, a similar phe¬ 
nomenology is expected for any depinning model as soon 
as the elastic propagator exhibits soft modes. As dis¬ 
cussed in [35] in the case of plastic yielding, this new 


sub-class of depinning model is expected to exhibit non¬ 
trivial scaling properties. More generally it is tempting to 
study in more details the ergodic behavior of such mod¬ 
els at finite temperature in relation with the Soft Glassy 
Rheology models m and with the recent observation of 
the strong effect of Eshelby events on relaxation processes 
in the liquid state [59] . 
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